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Abstract. Condensed matter exhibits a wide variety of exotic emergent phenomena 
such as the fractional quantum Hall effect and the low temperature cooperative 
behavior of highly frustrated magnets. I consider the classical Hamiltonian dynamics 
of spins of the latter phenomena using a method introduced by Dirac in the 1950s by 
assuming they are constrained to their lowest energy configurations as a simplifying 
measure. Focusing on the kagomc antifcrromagnet as an example, I find it is a gauge 
system with topological dynamics and non-locally connected edge states for certain 
open boundary conditions similar to doubled Chern-Simons electrodynamics expected 
of a Z2 spin liquid. These dynamics are also similar to electrons in the fractional 
quantum Hall effect. The classical theory presented here is a first step towards 
a controlled semi-classical description of the spin liquid phases of many pyrochlore 
and kagome antiferromagnets and towards a description of the low energy classical 
dynamics of the corresponding unconstrained Hcisenbcrg models. 
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1. Introduction 

There is now substantial experimental [U |2] and numerical [3J HI [2] evidence that a 
combination of frustration and low spin moments is the mechanism that produces novel 
quantum spin liquid phases in highly frustrated magnets (HFMs). It is not obvious, 
however, why this should be true. On the kagome lattice, large-iV spin models can 
exhibit order-by-disorder by quantum fluctuations at leading order in the semi-classical 
large S/N limit [6J. This induced order relieves the frustration of these "spins" and places 
them on similar grounds with unfrustrated ones at smaller S/N. So, given the evidence 
for disordered spin liquid phases, the order-by-disorder phenomena must be much more 
delicate for ordinary SU(2) spins. Such a conclusion is further backed up by heroic 
efforts to compute high order large-S expansions [7J [S]. To address the mechanism that 
produces these new spin liquid phases then we need a description of these materials 
that does not begin by relieving the frustration and/or relying on the dominance of 
order-by- disorder . 

To understand the connection between frustration understood at the classical level 
and the novel phases observed at low spin moments, one way to proceed would be to 
construct a semi-classical approximation that is still capable of describing non-ordered 
phases at smaller S. Such a description would be similar to the quantum melting of an 
antiferromagnetic phase due to order parameter fluctuations [9] only the dimension of the 
"order parameter" manifold, which are the classical ground states, would be larger, even 
growing with the size of the system in the highly frustrated case. At present, we know 
some features of this description. Similar to the distinctions between even, odd and half 
odd integer spinsfTlUfTTj in square lattice antiferromagnets, tunneling processes between 
classical ground states of kagome antiferromagnets leads to distinctions between integer 
and half odd integer spins[12j. However, it is unclear at present how these tunneling 
processes and other fluctuations produce a spin liquid phase at smaller spin moments. 

To make progress on such a semi-classical description of HFMs, here I study, using 
Dirac's "Generalized Hamiltonian mechanics" [T3 | [T ^ IT5 } ITE] . the "dynamics" of spins 
constrained to their ground state configurations. By counting the number of canonical 
coordinates needed to parametrize this surface, I show that it is a null surface for the 
Poisson bracket much like the light cone is a null surface in Minkowsky spacetime (See 
Ch. 2 of [IS]). In simpler terms, this means that the number of canonical coordinates 
of the constrained phase space is not equal to its dimension. Some coordinates are 
redundant non-canonical coordiantes. According to the Dirac method, this implies the 
constrained spin model has gauge dynamics. In particular, for the case of periodic 
boundary conditions, I find that the surface whose dimension is well known to grow 
with the system size[T7] has no canonical coordinates. For open boundary conditions, I 
find that this number depends on the existence of "dangling triangles" on the boundary 
of the kagome cluster and at best grows with its circumference. The constrained spin 
kagome model is therefore much like a "doubled" version of topological Chern- Simons 
electrodynamics characteristic of a Z2 spin liquid state [18]. However, it does not require 
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longer ranged interactions to construct as in the Levin and Wen models jl9J. This 
approach generalizes easily to other HFMs though the counting of canonical coordinates 
is likely different in other cases with, for example, the number likely growing with the 
volume of the system on the pyrochlore lattice. We will conclude with a discussion of 
the implication of these results for the ordinary unconstrained Heisenberg model and 
for more realistic models subject to various perturbations. 



2. Kagome ground state spin configurations 

In a highly frustrated magnet (HFM), spins are frustrated because they have many 
options to choose from and are unable to decide which is best. On the kagome lattice 
shown in Fig. [TJ the classical ground states of the nearest neighbor Heisenberg model 
prefer a vanishing total spin on each triangle [TT1 120] 

<f>i jk ,a = ^ + ^ + n a k = 0, ae {x, y, z}, (1) 

where a G {x, y, z} are the three components of the classical spin unit vector Oj on 
site i and the three sites i, j and k form any triangle on the lattice. This may happen 
in many materials including Herbersmithite, the Jarosite family, SrCrg^G^zOig 
and Na4lr308- Spins "suffering" this condition are highly frustrated for they have 
difficulty deciding between the continuously many arrangements that satisfy it. Such 
arrangements are described by the "spin origami" construction [2T1 122] of drawing 
spin vectors on a piece of paper and literally folding the paper to obtain new spin 
directions (see below). The resulting behavior of the spins is then collective and at 
finite temperatures, they enter a "cooperative" paramagnetic phase|23j. Furthermore, 
there is a wide class of other HFMs with similar constraints [24J such as the pyrochlore 
antiferromagnets where the analog of the spin origami construction leads to an effective 
Maxwell-like gauge description and dipolar spin correlations|25l |2T)] . Because of their 
novel low energy properties due to such constrained mechanical behavior, HFMs 
continue to be promising materials [27] to search for new phases of matter. 

In essence the origami sheet construction gives us an intuitive representation of the 
kagome ground state spin configurations. It is a duality transformation to a set of height 
vectors hj that live on the triangular lattice formed by the hexagons of the kagome 
lattice. These vectors are defined through the single spin shared by two neighboring 
hexagons and are directly analogous to displacement vectors of a two-dimensional solid 
membrane in three-dimensional space. Specifically, the mapping between height vectors 
and spin vectors is f2j = hi — hj where I J are the two hexagons that share the spin i. 
This naturally satisfies the constraint that Cli + Clj + Cl^ = on every triangle ijk of the 
kagome lattice provided a suitable convention for the signs of hi is made. Because f)j is 
a unit vector, however, we must always have \\hi — hj\\ = 1. A natural sign convention 
is to use the spin vectors in the q = configuration shown in Fig. [l'a) where the hi 
vector at the tip of a spin arrow enters with a positive sign and the hj vector at base of 
the arrow enters with a negative sign. Then setting the height of one of the hexagons to 
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Figure 1. The continuously many spin configurations of classical kagome 
antiferromagnet that have vanishing total spin on each solid triangle. The solid bonds 
are the kagome lattice, the dashed bonds the triangular lattice of hexagons, (a) The 
coplanar "q — 0" configuration viewed as arrows pointing to vertices of the dashed 
triangular lattice, (b) A "folding" of the spins along a dashed line of the triangular 
lattice viewed as a piece of paper, called a weather-vane mode[17j ETJ [22] • The folded 
spin directions also satisfies ([T]). All ground states modes arise from such folding of the 
spin paper, a construction called "spin origami" [21 j . They imply that low energy spin 
configurations evolve continuously and collectively in kagome antiferromagnets. (c) A 
visualization of the folded sheet in (b) using the solid analogy discussed in the main 
text. 



zero (say hi = (0,0,0)), all other heights can then be constructed recursively from the 
spin vectors such as h 2 = h 1 + Cl 1 , etc. The allowed height vectors are then obtained by 
literally folding the a piece of paper with the q = spin pattern drawn on it as in[T|b). 
The location of a point on this piece of paper is then hi and the arrows drawn on the 
paper are the spin vectors themselves. Fig. [I](c) is an alternative presentation of the 
folded paper in Figjljb) viewed as a two-dimensional triangular lattice sheet floating in 
a three dimensional space where each bond is exactly one unit in length (using Jmol, 
see 



http : //www . jmol . org). 



3. Degrees of freedom counting 

Given the complexity of the kagome ground state spin configurations, an important 
property to understand are the number of degrees of freedom. One measure of the 
number of degrees of freedom, discussed in Ref. [21] and frequently called Maxwell 
mode counting, is simply the number of free coordinates d = D — M where D is the 
total number of coordinates and M is the number of constraint functions such as those 
defined in 0. Naively, d = on the kagome lattice because D = 2N where N are the 
number of spins and M = 3N& where = 2N/3 are the number of triangles. Since 
we have already constructed several "folding" modes above proving d > 0, this naive 
argument must fail and so the constraints are not all independent (i.e. M < 3N&). 

Given this complexity, a better way to determine d is to look at the number of 
ways one can fold the spin origami sheet, i.e. its allowed folding patterns. The sheet of 
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Fig. WCb) clearly has two parallel folding modes so d = 2 for it. We cannot fold along 
any of the other modes without tearing the piece of paper and violating the constraint 
11^/ — hj\ | = 1. These two folds then define one folding pattern. However, if we were 
to flatten the paper, we could then fold along one of the other 6 lines. After this fold, 
we could no longer fold along the original two lines in Fig. [T^b), but would find that 
we could only fold along two other lines parallel to the one we chose without violating 
the constraint. Hence, this folding pattern has three allowed folds giving d = 3. The 
remaining three lines then make up a third folding pattern. So the flat sheet here, the 
coplanar q = state, is the intersection of three distinct folding patterns. These can 
alternatively be thought of as three smooth surfaces embedded in the unconstrained D 
dimensional coordinate space with dimensionalities d taken from {2, 3, 3} that meet at a 
point that defines the flat sheet. So d takes on different values depending on the origami 
sheet's folding pattern and the constrained space is a collection of intersecting surfaces 
defined by these folding patterns in a D dimensional space. 

When considering the dynamics of spins, however, the number of free coordinates d 
is not an important property. Instead it is the number of canonical coordinates N c that 
specify the dynamics. To understand this distinction better, lets demand, as we will 
throughout this paper, that each spin obey precessional dynamics. We can accomplish 
this by mapping the azimuthal ip and polar 6 coordinates of the spin unit vectors onto 
a position q and momentum p variable, £l(<p, 0) — > Cl(q,p) such that they obey the usual 
angular momentum relations {Q x ,Q y } = Q z /S where S is the spin length or quantum 
number and {, } is the usual Poisson bracket. One choice is q = ip, p = S cos8 = Sfl z . 
In an unconstrained spin system, we then see that N c = D = 2N. However, if we 
were to impose an odd number of constraints M we clearly cannot have d = N c for 
N c must be even. More generally, Dirac found by using Lagrange multipliers to impose 
the constraints that iV c <= d and that when they are unequal, the extra Nq = d — N c 
coordinates are redundant as far as time evolution is concerned and naturally thought 
of as gauge degrees of freedom. 



3.1. Triangle and bow-tie models 

To place the above discussion in a simple context relevant to the dynamics of kagome 
antiferromagnets, let us find N c for the simpler triangle and bow tie systems, shown 
in Figs. [2^i and [3^,, systems with only one folding pattern, before turning to the full 
kagome lattice system. Additional calculations of this sort are provided in Appendix 



|A| where N c is calculated for several well known gauge systems, to place these sort of 
calculations in a more familiar context. 

The first step to determining N c is to find d. The unconstrained phase space of the 
triangle system, that describes the configurations of the three spins ^(g 1 ,^), fi(g 2 ,p 2 ) 
and Cl(q 3 ,p 3 ), clearly has D = 6. Imposing the M = 3 constraints (f>i23, x = (pi23, y = 
^i23,z = then tells us d = 3 since they are independent constraints. Because d is odd, 
there is necessarily some ambiguity in identifying canonical degrees of freedom and so 
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Figure 2. Dynamics of the triangle system, (a) The three constrained spins of the 
single triangle that lie on the plane defined by the dashed triangle origami sheet. All 
spin configurations are then global rotations of the form Si = H-S®, that can be viewed 
as the dashed triangle floating in three dimensional space, (b) Two trajectories related 
by a deformation plotted in the well known solid ball that parameterizes R € 50(3). 
For each choice of h\23{t), a different trajectory will result from solving the equations 
of motion starting from a given initial condition labeled by i. If two choices differ 
by a smooth function, then we can view one trajectory as a deformation of the other 
much like a coffee mug and a doughnut can be deformed into each other. Since all 
trajectories are deformable into each other through these different choices of h\2z(t), 
this system has only one possible trajectory. 

this system must have the gauge dynamics discussed above, its phase space must involve 
at least one gauge coordinate. 

The next step to determining N c it to then work in the unconstrained phase space 
by introducing Lagrange multipliers and use them to impose the constraints. We do so 
by extending the Hamiltonian to He = H—^2 ia hf 2 34>i23,a where the Lagrange multipliers 
take the form of magnetic fields h\23- We then choose h\2z so that 0i23, a = so that if we 
obey the constraints at time t = 0, we will do so for all future times. For a Heisenberg 
model H = JJ2(ij) ^ * I obtain 

<fil23,x = {<Pl23,x, H}-^2 hl23{<f>123,x, 4>123,a} = (2) 

a 

independent of hus provided initially 4>i23, a = . Similarly 0i23, y = ^123,2 = 
is independent of /ii 23 . So there are N L = 3 Lagrange multiplier functions ^123 
undetermined by imposing 0i23, a = at t = 0. For each choice of ft. 123, we will then get 
a different time evolution or trajectory for Cli that nevertheless obeys the constraints 

0123,a = 0! 

The final step then is to make sense of the above result. Given a choice of 
the arbitrary Lagrange multipliers functions ^123 (i), viewed as a "choice of gauge", a 
trajectory in this phase space follows from the equations of motion such as the cartoon 
picture of [2^. If we smoothly change hi2s(t) we would obtain a new "gauge equivalent" 
trajectory that can be viewed as a deformation of the first as in [2j3. The study of the 
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Figure 3. Dynamics of the constrained bow tie spin system, (a) Origami construction 
of the spins in the bow tie system viewed as spins drawn on the dashed triangles. All 
spin configurations are spanned by a choice of f^, the folding angle 9 and an angle <p 
defined by global spin rotations around (b) Trajectories in the four dimensional 
constrained phase space of the bow tie system projected onto the torus formed by the 
folding angle 9 and the angle ip starting from a reference spin configuration. Under time 
evolution, only the two (gauge) coordinates 9 and <p depend on the different choices of 
hi23(t) and /i34s(i). Any two trajectories with the same &3(t), such as the two shown, 
are deformable into each other by changing hi2z{t) and h^(t). 



topology of a set works in a similar way. To see that a coffee mug and a doughnut are 
topologically equivalent we deform the mug into the doughnut. Here, utilizing different 
choices of /1123(f) we can deform trajectories into other trajectories. Since there are 
Nl = 3 of these Lagrange multiplier functions the space of trajectories we obtain is 
necessarily three dimensional and spanned by Nq = Nl = 3 gauge coordinates. Hence 
we obtain Dirac's formula for the counting of canonical degrees of freedom: iV c = d — Ni 
and obtain N c = for the triangle system that reflects the absence of any freedom to 
have different trajectories that are not related by different gauge choices of /2.123(f)- 

Other examples of systems with similar dynamics include a charged particle in 
a very large magnetic field with a fixed angular momentum [25] and "Chern-Simons" 
electrodynamics in two spatial dimensions [29] (both discussed in Appendix A). The 
latter example is perhaps the most well known for it is related to both Einstein gravity 
in two spatial dimensions [30] and the fractional quantum hall effect (FQHE) of a two 
dimensional electron gas in a very large magnetic field[3lj, [32J [33] . Though Chern- 
Simons electrodynamics has no canonical coordinates like the triangle model, it does 
have a form of non-trivial dynamics. 

In the bow-tie system of two triangles (M = 6) and five spins (D = 10, d = 4), 
shown in Fig. [3j the counting of arbitrary Lagrange multipliers demonstrates that it 
has two canonical and two gauge coordinates. To see this, we set 0i23, a = 0345,a = 
and solve for /1123(f) and /1345(f). Since any choice with fields parallel to Cl 3 , i.e. 
ft- 123 (t) = ft, 123(^)^3 and /1345(f) = ^345(^^3, satisfies these equations, there are two 
arbitrary Lagrange multiplier functions giving Nl = 2 and hence N c = 2. Also, by 
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inspection, we can identify the stated two canonical degrees of freedom. These are the 
two coordinates needed to specify Cl 3 because they evolve in time independent of the 
choice of lengths hi23 arid /1345. After fixing Q 3 , two additional angles are required to 
specify a configuration of the bow tie system: the "folding angle" 9 shown in Fig. [3^i 
and an angle ip that determines the global rotation about Cl 3 of all spins from a chosen 
reference configuration. As in the triangle model, we can again make a choice of hi23(t) 
and ^1345 (t) and compute a trajectory from the equations of motion. Here we find, by 
smoothly changing these two functions, that a two dimensional space of trajectories can 
be deformed into each other all with the same value of f2 3 . Namely, Cl 3 is gauge invariant 
and its two canonical coordinates span the space of trajectories not deformable into each 
other through changes in h^if) and /1345 (i). Hence, the dynamics is completely specified 
by the behavior of Cl 3 and this five spin constrained system reduces to the description 
of a single spin. 

4. Full kagome lattice system 

Consider now the full kagome lattice system. To count and study its canonical degrees of 
freedom, I have implemented a computational scheme based on the local properties near 
one spin configuration in the constrained phase space as discussed below. A selection of 
the counting results are then presented in table [1} For periodic boundary conditions, I 
find N c = for any sized system, including the smallest system size that is equivalent to 
and in agreement with the single triangle system discussed above. However, given that 
the constrained space of the Full kagome lattice is a collection of intersecting surfaces 
as discussed in [3j it is not obvious whether all trajectories can be deformed into each 
other like in the triangle model. Hence the dynamics of the constrained spins in the 
kagome antiferromagnet with periodic boundary conditions is composed of one or more 
discrete sets of topologically equivilant trajectores. 

4-1. Counting Algorithm 

To study the degrees of freedom of the full kagome lattice model, let us impose ((!]) 
on every triangle (ijk) of the lattice. To determine the number of arbitrary Lagrange 
multiplier functions, we extend the Hamiltonian to He = H — Yliijk) ^ijk ' <<% fc anc ^ sorve 
</>ijfc,a — so that if at time t = 4>ijk, a = 0, it remains zero for all time. Here, <pijk, a 
is given by the usual Poisson bracket relation for the time evolution of any phase space 
observable: 

4>ijk,a = {4>ijk,a, He} = {4>ijk,a, H} — hl m n t b{4>ijk,a, 4>lmn,b} , (3) 

(lmn),b 

where (ijk) and (Imn) denotes a triangle, a and b range through {x, y, z} and we have 
assumed 4>ijk,a = at £ = 0. For the Heisenberg model 

iJ=JS 2 X>-% = ^ + const. (4) 

(ij) (ijk) 
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we have {<pijk, a , H} = for initial conditions with 4>ijh,a = 0. Viewing {4>ijk,ai (fiimn,b} as 
a square antisymmetric matrix C ai p, with a <H- ijk, a and (5 -H- Imn, b, we are then left 
with the following eigenvalue problem 

p 

If C aj /3 is invertible, then this equation has a unique solution of hp = and no arbitrary 
Lagrange multipliers (Nl = 0). However, if C a ^ is not invertible, if it has zero 
eigenvalues, then it has many solutions spanned by the null space of C a $. Suppose 
we find one solution h a = haK By knowing an eigenvector X a of with zero 
eigenvalue, ha = ha + x(t)X a is also a solution for any x(t). This function x(t) is 
then an arbitrary Lagrange multiplier. So to determine N c we need only determine the 
number of zero eigenvalues of C a ,p, i.e. the dimension of its null space, for this equals 
the number of arbitrary Lagrange multipliers Nl and N c = D — M — Nl with M the 
number of constraint functions <pijk, a and D = 2N as before. 

The above describes in a nutshell the algorithm for counting the number of canonical 
degrees of freedom. However, there are two subtleties that need to be considered before 
applying it. The first is to construct C a ,/3 out of an independent set of constraint 
functions. Given the form of He, if some 4>ijk,a could be written as a linear combination 
of the others, it is redundant and should be removed from the set of constraints. We 
will find this happens frequently. The second subtlety is to avoid computing C a ^ at a 
coplanar spin configuration or other point where several distinct folding patterns meet 
as discussed in |3j This is easily achieved by choosing a generic spin configuration. 

To find a linearly independent set of constraint functions, lets define a set of vectors, 
one for each constraint, by taking the phase space gradient of each constraint function, 
v = (d<p a /dq l , d(p a /dpi, . . .) evaluated at some point (q l ,pi,...) in the constrained 
phase space. We then need only chose a linearly independent set of vectors w 5 , 
a = 1 . . . M using the Gram-Schmidt procedure to obtain M. The constraint matrix of 
the independent set is then readily evaluated through, 

= W &,1 W P,2 - w a,2W- p i + Wa&Dfa ~ ^5,4 W /3,3 + ( 6 ) 

The number of zero eigenvalues of this matrix is then Nl as discussed above. 

A python script implementing this counting algorithm, available online, is discussed 



in more detail in Appendix B The results for various system sizes and boundary 
conditions obtained at multiple generic points on smooth portions of the constrained 
phase space are presented in [Tj Remarkably, the results for iV c were always the same 
independent of which folding pattern the generic point belonged. So unlike d, the 
dimension of a given folding pattern that changes from one pattern/surface to another, 
the number of canonical coordinates is always the same. 

For open boundary conditions, I find N c > and evidence for unusual edge states. 
These are similar to the edge-states responsible for the vanishing longitudinal and 
quantized hall resistance that define the FQHEJ31]. As shown in table [TJ N c is equal to 
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State 


Lattice size 


Boundary Conditions 


D 


M 


Nr 


± * C 


Near a/3 x a/3 


3x3 


periodic 


54 


49 


5 





Near g = 


3x3 


periodic 


54 


49 


5 





Near g = 


4x4 


periodic 


96 


90 


6 





Near q = 


5x5 


periodic 


150 


143 


7 





Near a/3 x a/3 


3x3 


open 


70 


53 


11 


6 


Near g = 


3x3 


open 


70 


54 


10 


6 


Near g = 


4x4 


open 


118 


96 


14 


8 


Near q = 


5x5 


open 


178 


150 


18 


10 


Near g = 


3x3 


cylindrical — no A 


60 


51 


7 


2 


Near g = 


4x4 


cylindrical — no A 


104 


92 


10 


2 


Near g = 


5x5 


cylindrical — no A 


160 


145 


13 


2 


Near q = 


3x3 


cylindrical — 6 A 


66 


54 


6 


6 


Near g = 


4x4 


cylindrical — 8 A 


112 


96 


8 


8 


Near g = 


5x5 


cylindrical — 10 A 


170 


150 


10 


10 



Table 1. Counting of degrees of freedom near two spin configurations and for various 
lattice sizes and boundary conditions. Here D is twice the number of spins and 
the number of unconstrained degrees of freedom, M is the number of independent 
constraint functions and Nl the number of arbitrary Lagrange multipliers. Note: 
"cylindrical — no A" means cylindrical boundary conditions without dangling triangles 
and "cylindrical — 6 A" means with six dangling triangles. 

the number of "dangling triangles" (see Fig. 13(a)) unless there are none in which case 
N c = 2. The bow tie system, with two dangling triangles and N c = 2, is the simplest 
example of this result. Identifying the canonical degrees of freedom with the dangling 
triangles leads to a remarkable conclusion. Because two degrees of freedom, a position 
and a momentum variable, are needed to describe a local mechanical object and there 
is only one canonical degree of freedom per such triangle, they must be non-locally 
connected; this result is in agreement with the direct study of the trajectories discussed 
in the next subsection. 

4-2. Orbits in the constrained phase space and edges modes 

One way to understand and identify the physical degrees of freedom is to construct 
the "reduced" phase space of the constrained phase space. That is, if we identify all 
coordinates that are redundant we can then ignore them and focus on the remaining 
canonical coordinates. The reduced phase space solves the gauge-redundancy problem 
explicitly. This solution is particularly useful for our purpose of identifying canonical 
coordinates numerically. 

A nice way to understand the reduced phase space is to construct a basis of tangent 
vectors at a point (q l ,Pi, . . .) in the constrained phase space and find the subset of these 




Figure 4. Canonical degrees of freedom of the constrained phase space, (a) An 
example of boundary conditions with six dangling triangles highlighted at the top and 
bottom and N c = 6. Left and right sides may or may not be connected to form a 
cylinder, (b)-(d) Examples of orbits visualized using the origami sheet construction 
with Jmol (http://www.jmol.org). The white bonds correspond to the external 
dangling triangle spins, the blue bonds the third spin on each dangling triangle and 
the green bonds the bulk spins. Starting from (b), the initial spin configuration, (c) 
is a mode that "folds" the left side and (d) is the conjugate mode to (c) with opposite 
behavior on the top and bottom blue edge bonds. Notice how more than one dangling 
triangle move for both modes, an observation that is generally true of all orbits that 
changing the canonical degrees of freedom. 



associated with redundant coordinates. The remaining vectors point along canonical 
coordinate axes. It turns out that if X a is an eigenvector of C a ,/3 with zero eigenvalue 
then the phase space vector 

G = Vl (— -— — -— ^ (7) 

is tangent to the constrained phase space and points along a gauge coordinate axis 
known as a "gauge orbit" . The dot product of this vector with phase space gradient of 



a constraint function v a (defined previously in subsection 4.1) vanishes for 



But the vectors v a , as gradients of the constraint functions, are a complete set of normal 
vectors to the constrained phase space surface. So, from the eigenvectors of C a ^ with 
zero eigenvalues, we can immediately construct tangent vectors. A complete basis of 
tangent vectors can then be constructed by finding all vectors Pj perpendicular to v Q and 
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adding as many of them to the set {G} as we can to obtain a set {Gi, G2, . . . , Pi, P2, • • •} 
that forms a basis of the tangent space at the point . . .)■ 

As stated, the vector G points along a gauge orbit; the other Pj vectors then point 
along "physical" directions. To see this, consider the equation of motion for any phase 
space observable 

f = {f,H}-J2hMAa} (9) 



4.1 



we obtain 



Letting h a = ha + x(t)X a , as discussed in section 

f = {f,H}-J2 h${f, «} ~ x{t) X M, «} (10) 

a a 

= {/, H}-J2 h${f, 4a) - x{t)F ■ G (11) 

a 

= {f,H'}-x(t)F-G (12) 

where F = (df /dq 1 , df/dpi, . . .) and H' = H — ha (f> a - So G is directly associated 
with the arbitrary Lagrange multiplier function x(t) and any observable / with F 
perpendicular to G evolves in time independent of this x(t). Now suppose / was the 
phase space coordinate q 1 and we evolve it in time starting from a given initial conditions 
g x (0) under two different choices of x = x\ and x = x 2 - At an infinitesimal time t = e 
later, the two new points g x (t = e; xi) and g 1 (t = e\ x 2 ) are related by 

q 1 (e\ x x ) - q 1 (e\ x 2 ) = e(x x - x 2 ){l, 0, . . .) • G (13) 

where (1, 0, . . .) • G is just the component of G along the q 1 direction in phase space. So 
performing this same calculation for the other coordinates, we see that nearby physically 
equivalent points in phase space are connected by the vector G and path from one of 
these points to the other is an unphysical gauge orbit. The remaining tangent vectors 
Pi, P2, • • • that are linearly independent from G then point along directions that cannot 
be related by a change in x(t). 

Lets now apply this vector analysis to show that trajectories unrelated by a change 
in x(t) necessarily involve the motion of spins on more than one dangling triangle. To 
this end, we need to show that all linear combinations of Pi, P2 etc., move spins on more 
than one dangling triangle. For a given dangling triangle, we then need only define a set 
of vectors F = (0, . . . , 1, 0, . . .) corresponding to the q % or pi coordinate of a bulk spin 
or a spin on our chosen dangling triangle and show that the set {F} + {P} is linearly 
independent. For the case shown in Fig. |4| the linear independence of this set was 
found to be true for each dangling triangle. Hence, the canonical modes, as expected 
from the inability to assign both a position and momentum variable to the spins on 
a given dangling triangle, are non-local and necessarily involve the motion of spins on 
more than one dangling triangle. 

It is useful to use the spin origami sheet construction to visualize the orbits 
defined by a tangent vector Pj. Using this vector, we can construct a small orbit 
by y(t) = y(0) + tP, < t < e. Here y = {q l ,P\, ■ ■ .) is a point in the unconstrained 
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phase space. This orbit essentially moves a small distance away from y(0) and shows 
how physical modes evolve in time in a given gauge. To see how this affects the spins, 
it is then straightforward to map y to the set of spin vectors f2j. It is not easy to 
understand the collective behavior by directly observing individual spin vectors so lets 
pass to the spin origami construction for visualization purposes. The result is presented 
in Fig. [4] b-d and movies of them are available online in the supplementary materials. 
The motion of multiple dangling triangles (that have white or blue bonds) is apparent 
in both the movie and in Fig. |4](b)-(d). 

5. Discussion 

In summary using the Dirac approach, I found that the number of canonical coordinates 
N c , unlike the d — D — M coordinates given by Maxwell counting, is the same on all the 
intersecting surfaces characterized by folding patterns that make up the constrained 
phase space. In addition, this N c vanishes on any cluster with periodic boundary 
conditions and grows with the number of dangling triangles on clusters with open 
boundary conditions. The remarkable simplicity of this result allows us to make 
connections both with known results for the quantum Heisenberg model and for systems 
whose Hamiltonian is perturbed away from the nearest neighbor model. 

The results for N c here are very similar to doubled Chern-simons electrodynamics 
expected to describe Z 2 spin liquids[T8]. The Z 2 spin liquid picture recently gained 
greater acceptance through DMRG calculations showing a gap to all bulk excitations [5] 
and entanglement entropy calculations consistent with a topological Z 2 spin liquid [35, 
136] . This suggests that even spin 1/2 spins may know about the classical constrained 
phase space that forms the focus of this paper! 

The canonical degree of freedom counting also has implications for realistic 
Hamiltonians of the form 

H = J ^2 Si ■ Sj + K (perturbations) (14) 

where Si are quantum spin operators, the first term is the nearest neighbor exchange and 
the second, characterized by an energy scale K that is much smaller than J, represents 
all other perturbations, including further neighbor exchange, Dzyaloshinskii-Moriya 
interactions, ring exchange and impurities. In the limit J — >■ 00, the perturbations 
would induce a Hamiltonian for the constrained phase space with energy scale K. Likely 
this induced Hamiltonian would lead to so-called "secondary constraints" and freeze the 
spins into a single pattern, their ground state. For finite J, however, where both gauge 
and canonical modes are "physical modes", the perturbing Hamiltonian would lift the 
degeneracy of these zero modes giving them dispersion. However, because these modes 
transform differently under gauge transformations in the constrained model, they should 
remain as two distinct types of modes. The low energy sector is then spanned by two 
kinds of zero-modes (the gauge and canonical modes of the constrained phase space) 
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and a third mode involving fluctuations outside of the constrained space such as the 
"monopoles" of Ref. [26] . 
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Appendix A. Physical degrees of freedom counting for several classical 
gauge systems 

Maxwell Electrodynamics 

It is well known that electromagnetic waves are transverse with only two polarizations. 
Viewed as a Hamiltonian mechanical system, each polarization mode should consist 
of a position and momentum variable for each point in space. This means that 
electrodynamics has four canonical degrees of freedom per point in space. However, 
in the Lagrangian view of electrodyamics (without source terms), the action 

S Maxwell = J d^C MaxwM = -\J d^F^ (A.l) 

with F^y = d^A u — d u A^ is described in terms of four fields = (At, A), the scalar 
and vector potentials. For each "position" field A^ there should be a "velocity" field 
dtA^ making a total of 8 degrees of freedom per point in space. From this perspective, 
it is not obvious that there are two polarizations of light. Lets show this by passing 
to the Hamiltonian formalism. A discussion along these lines is also available in Ref. 
[15] but here I present it in the language of the main body of this paper to facility an 
understanding of the calculations leading to its main result that N c = in the bulk of 
the kagome constrained spin model. 

To construct a Hamiltonian, we first need to find the momentum fields 7r M = 
SS Maxwell /Sd t Afj,. They are 

^ = d n A t _ d t A „ ( A _ 2 ) 

The time component of this equation, 7r* = <pi = is a phase space constraint, for it 
does not relate the momentum variable 7r* to any velocity variables dtA^. The spatial 
components 7r = — VA — d t A = E, tells us that the electric field components are the 
other momentum variables. This means that Gauss's law V ■ E = <p2 = is also a phase 
space constraint for it is only a relationship between the momentum densities jf. So, 
the constrained phase space is parametrized by the eight variables 7r M and A^ subject 
to the two constraints <f>\ = and 02 = 0. 

To find the number of physical degrees of freedom per point in space, N c , the 
approach used in the main text is to work in the unconstrained phase space using 
Lagrange multipliers. Then N c per point in space is given by N c — D — M — Nl where 
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D = 8 is the number of unconstrained degrees of freedom, M = 2 is the number of 
constraint functions and Nl the number of arbitrary Lagrange multipliers not fixed 
by the requirement to remain in the constrained phase space. This requirement is 
constructed by starting from the Hamiltonian, 

U = J d 3 r[7r^ - £] = iy d 3 r(7? 2 + (V x A) 2 ), (A.3) 

extending to He = H + f d 3 r [Ai</>i + A202] and choosing Ai and A2 so that 

<Mr) = {Mf)M + / <Pr%(?){<h(f),<h(?)} = 0> ( A - 4 ) 
2 (r) = {0 2 (r),^} + / d 3 r'X 1 ^){<f )2 (r),Mn} = 0- (A.5) 



All Poisson brackets in these expressions vanish so both Ai(r) and A2(r) are arbitrary 
and iVx, = 2. Hence A^ c = 4 per point in space and as discussed above these four degrees 
of freedom correspond to the two polarizations of light. 

Abelian Chern-Simon's theory 

Abelian Chern-Simon's theory is formally very similar to electrodynamics in two spatial 
and one time dimension. Maxwell electrodynamics in this two-dimensional case would 
still have a Lagrangian density —\F ia ,F ttv , just \i G {t,x,y} does not include the z- 
direction. However, with only two spatial dimensions, the action 

k f 

Schern-Simons = ^ 0? r dtE^ X A v A X (A.6) 

is also allowed. Notice it does not involve the "metric tensor" g^ u but instead only 
the antisymmetric tensor e fluX . This is a hint that space-time distance, defined by 
s 2 = giiv^x" = t 2 — x 2 — y 2 — z 2 , may not be important in computing this action but that 
it may only depend on the topology of the space it is integrating over. It is therefore an 
interesting alternative to electrodynamics that can arise in a lower dimensional setting. 

To count the physical degrees of freedom of Chern-Simons electrodynamics, let us 
again pass to the Hamiltonian formalism using the language of the main text. See also 
Ref. [37]. The unconstrained phase space has D = 6 degrees of freedom per point in 
space with the three momentum variables 

77* = 0, ir x = — A v , n y = -—A x . (A.7) 
2tt y ' 2tt v ' 

Because none of these three equations involve the velocity variables dtA^ they are all 
phase space constraints. Let us define these constraints through the functions 

k k 
(f> t = n\ (f) x = tt x - —A y , (f) y = ir y + —A x (A.8) 
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So we have three constraints 0^ = from the definition of the momentum variables. To 
see if there are any other constraints, we need to look at the equations of motion. They 
turn out to be the three equations 



2n ^ 



2tt 



(d„A v - d v A») = 



(A.9) 



so that there are no electric or magnetic fields. One of these equations, 04 = 
j- [d x A y — d y A x ) = 0, does not involve any time evolution. It is therefore also a 
constraint. 

In total then, we have M = 4 constraints, M and 04. Because there are no electric 
and magnetic fields, the Hamiltonian H = vanishes. Introducing Lagrange multipliers 
through 

H E = H + J d 2 r [A^ + A 4 4 ] , (A. 10) 

we need to solve 



r),0 b (r<)}AV 



J2 I d 2 r'C ab (ry)X b (n = (A.ll 



for X a (f) where b G {t, x,y, 4}. One solution is just A a (r) = 0. All solutions define 
the null space of the constraint "matrix" C a b(f, r 1 ) and they are spanned by the 
eigenfunctions of C a fe(f*, r*) with eigenvalue 0. Computing the matrix explicitly, I find 



C(r, r) 



(o 








\ 







7T 


2k U V 







i o 


_Jlq> 


5 




-JLff JLff 
2k U V 2-k u x 





/ 



(A.12) 



Fourier transforming and diagonalizing the resulting 4x4 matrix, we find two 
eigenvalues. Transforming back to real space, we discover the null space is spanned 
by the orthogonal eigenfunctions 



A*(f) = ^i(f), X x = X y = X 4 







(A.13) 



and 



A*(f) = 0, X x (f) = ^ 2 (r), \y(r)= l -d y ^ 2 {r), A 4 (f) = ^ 2 (r), (A.14) 
where T/>i(f) and ^(r) are arbitrary functions. These two solutions can be verified by 



direct substitution into (A.12). Given the two arbitrary Lagrange multiplier functions 
■01 (r) and ^(^0 we have -^Vl — 2 per point in space. Hence, the number of canonical 
degrees of freedom per point in space is N c = D — M — Nl = 6 — 4 — 2 = 0. There are 
no canonical degrees of freedom in Chern-Simons electrodynamics. 
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Simple analog of Chern-Simons theory in classical mechanics 

In Ref. [2H], "Topological" (Chern-Simons) quantum mechanics, by Dunne, Jackiw 
and Trugenberger, a simple mechanics model was introduced through an analogy with 
Chern-Simons gauge theory. Here, using the language of this paper, I will show how 
similar in structure this model is to the kagome constrained spin model discussed in the 
main text. 

Ref. [2B] considers the mechanics of charged particles subject to an external 
magnetic field in two dimensions (inspired by the fractional quantum hall effect). In the 
circular gauge, the Lagrangian of such particles is 

L = ■ q + eA(q) ■ q = y {q 2 x + q 2 y ) + ^ (q x q y - q y q x ) (A. 15) 

In the fractional quantum hall effect, the important limit of the above model is that 
of very large magnetic fields B — » oo. As such, the authors consider the limit m — > 0. 
After taking this limit, the Lagrangian is linear in q and the conjugate momenta are 

e B eB . . 

Px = Py = ~yQx (A. 16) 

So there is no relationship between the momenta and the velocities. In passing to the 
Hamiltonian formalism, these M = 2 equations are the constraints 

/ e B , eB . 

(Px = Px + ( Py=Py- -J-Qx- ( A -!7) 

in the D = 4 dimensional unconstrained phase space and the Hamiltonian is H = 
Px ■ Qx + Py " Qy — L = 0. We can work in this phase space by introducing Lagrange 
multipliers through H E = H + u x (p x + u y (p y and imposing 

<j>i = {(pi, 4>j}uj = CijUj = (A. 18) 

where summation over j is implied. The matrix C = eB ( _°i o ) is invertible so there 
is one unique solution to this equation m = and there are no arbitrary Lagrange 
multipliers. This system is not a gauge system. The number of degrees of freedom, all 
of which are canonical, is N c = D — M — Nl = 4 — 2 — = 2. 

The focus of the paper, however, is not directly on this model but on the extension 
of this model to include an additional constraint, that of fixed angular momentum. This 
alternative model is described by the Lagrangian 

eB 

L = —(q x (q y - aq x ) - q y (q x + aq y )) + va (A. 19) 

where v is a parameter and a is a Lagrange multiplier that enforces a third constraint 

dL eB 9 , . s 



so that this system has M = 3. This fixes the angular momentum J because after 

eB 
2 



imposing cj) x = and <p y = 0, J = q x p y — q y p x = ^f-q 2 . Following the previous discussion 
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on fixing the Lagrange multipliers, we construct the constraint matrix C a b = {4>a,4>b}, 
a = {x, y, J} and obtain 



C 



/ eB -eBq x \ 

—eB —eBq y 
\eBq x eBq y J 



(A.21) 



This matrix has eigenvalues 0, ±i\/ B 2 + 2uB and the length of the eigenvector that 
corresponds to the zero eigenvalue 



u , 



Qx, u y = q y , uj = 1 (A.22) 



In [2] 
In [3] 
In [4] 



is arbitrary. So this model, with the additional constraint on angular momentum, has 
one arbitrary Lagrange multiplier and N c = 4 — 3 — 1 = physical degrees of freedom. 
It is very similar in structure to Chern-Simons electrodynamics, but in the context of 
the mechanics of a particle rather than a field. 

Appendix B. Notes on included python scripts 

There are two python scripts available online. To use them, you need scientific python or 
pylab. The script "state. py" holds the information about the spin configurations and 
lattices and the script "constraints .py" implements the counting scheme discussed in 
section 4.1. A typical use of these scripts would be 

In [1] : import state 

import constraints as con 
y_cs, T = state. Tl() 
state . check(y_cs ,T) 
Sum of deviations: 4 . 4408920985e-16 
In [5] : con. report (y_cs,T) 
No unconstrained degrees of freedom D: 6 
Total number of constraint functions: 3 
No of independent constraints M: 3 
Dimension d=D-M of the constrained phase spcae: 3 
Number of arbitrary Lagrange multipliers N_L: 3 
No of physical degrees of freedom: 

Here y_cs is the list of q's and p's associated with a given spin configuration 
{Q(q 1 ,pi),Q(q 2 ,p2), . . .} and T is an array that lists the sites in each triangle. In this 
case, T is just array ([ [1 ,2,3] ] ) because this example is the single triangle system 
created by the function state. Tl(). The call to state . check (y_cs ,T) is to make 
sure the constraints are obeyed by the spin configuration y_cs. To check the bow tie 
example, use the function state .T2() (in this case, T is array ( [[1,2,3] , [3,4,5]])). 
To study the full kagome lattice, six situations are included. Their boundary 



conditions for the 3x3 unit cell case are explicitly presented in Fig. Bl 
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Figure Bl. Explicit representation of the four boundary conditions used in the 
code for the 3x3 unit cell case, (a) Periodic boundary conditions, (b) Open 
boundary conditions, (c) cylindrical boundary conditions with no dangling triangles, 
(d) cylindrical boundary conditions with dangling triangles. 

• state . KqOpbc (Nl , N2 , dtheta) : This function looks at the N\ x N 2 unit cell lattice 
with periodic boundary conditions where each weather-vane mode is rotated away 
from the u q = 0" configuration by dtheta, 2dtheta, 3dtheta, etc. If dtheta is not 
specified, then it is set to the small irrational number \/2/100. 

• state . KqOobc (Nl , N2 , dtheta) : This function looks at the N\ x N 2 unit cell lattice 
with open boundary conditions near the u q = 0" spin configuration. 

• state . KqOcbcl (Nl , N2 , dtheta) : This function looks at the Ni x N 2 unit cell lattice 
with cylindrical boundary conditions that leave dangling triangles intact by wrap 
along the left and right side edges. It is also near the u q = 0" spin configuration. 

• state . Kq0cbc2 (Nl , N2 , dtheta) : This function looks at the N± x N 2 unit cell lattice 
with cylindrical boundary conditions that leave no dangling triangles by wrapping 
along the top and bottom leaving the left and right edges free. It is also near the 
"g = 0" spin configuration. 

• state . Kr3xr3pbc (dtheta) : This function looks at the 3x3 unit cell lattice with 
periodic boundary conditions near the "V3 x v^" spin configuration. 
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• state. Kq0pbc(Nl,N2,dtheta): This function looks at the 3x3 unit cell lattice 
with open boundary conditions near the "v^ x ^3" spin configuration. 

These routines in conjunction with con . report (y_cs , T) produced the results presented 
in Table [TJ 
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